PCX 



WORLD INmLECrUAL PROPERTY ORGANIZATION 
International Bureau 




INTERNATIONAL APPUCATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International PirtaitClMS^ 
G06F157/D0 



Al 



(U) InternatkHial Publkatoi Nnmbo^ WO 97/D7475 

(43) fiotematkmal Pablicatifm Date: 27 Febniaiy 1997 (27^97) 



(21) bttcniatimial ApplkatkA Number: Pa7US9S/10363 
02) biternaticmal FQiiig Date: 15 Aiigust 1995 (15j08^ 



(71) Appfieant: THE TRUSTEES OF COLUMBIA UNIVERSITY 

INTHEOTYOTNEW Y(»K[US/US}; n6th Street and 
Broadway. New Yoik. NY 10027-6699 (US). 

(72) Inventors: TRAUB. Joseph, F4 Coluinbia UnWemly, Com- 

piler Science TkpL, New Ycak, NY 10Q27-6699 (US). 
PASKOV. Spassimir, Cohmilna Univeisity, Conqmter Sci- 
ence DepL, New Ycnk. NY 10027-6699 (US). 

(74) Agents: TANG, Hemy, Y., S. ct aL; Bnimbau^ (haves, 
Donohue & Raymond, 3D Rock^dler Plaza, New Y«k, 
NY 101 12 (US). 



(81) Dedgnated States: CA, JP. KR, SG, European patent (AT, BE, 
CH, DE, DK. ES, FR, GB, OR, IE, IT, LU, MC, NL, PT, 
SE). 



Publi^ied 

With intemationai search report 



(54) Tide: ESTIMATION METHOD AND SYSTEM FOR FIN ANOAL SECURITIES TRADING 
(57) Abstract 



In siting tbe initial oHiering pnce of a financial instrument for purposes of 
securities trading, or in later revaluation as economic factms, sudi as interest rates, 
may change, an estimate of the vahie of the instrument may be represented as a multi- 
dimensional iutcgraL For evaluation of tbe inti^ral, numerical integration is preferred 
with tbe integiand being sampled at detcrministtc pcunts bavti^ a low-discicpancy 
property. The tcctoiique produces a p pro xim ate values at significant computational 
savings and wift greater idiabiHty as co m pared widi the Monte-Carlo technique (11). 
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Estimation Method And System For 
F;<^?"gj,?'' securities Trading 

The United States Government has certain rights 
to this invention pursuant to grants CC»- 91 -14042 cuid 
IRI-92-12597 awcurded by the National Science 
Foundation, and to grant AFOSR-91-0347 awarded by the 
5 U.S. Air Force. 

Background of the Invention 

The invention relates to financial securities 
trading such as, e.g., trading in stocks, bonds and 

10 financial derivative instruments , including futures, 
options and collateraJ.ized mortgage obligations. 

In financial securities trading, which includes 
the initial offer for sale, the voULue of a security may 
be estimated, e.g., based on eaqpected future cash flow. 

15 Such cash flow may depend on variable interest rates, 

for exanple, and these and other relevant variables may 
be viewed as stochastic variables. 

It is well known that the value of a financial 
security v^ch depends on stochastic varicdDles Ccin be 

20 estimated in terms of a imilti -dimensional integral. 
The dimen sion of such an integral tnay be very high. 

In collateralized mortgage obligations (CMO) , for 
example, instrtnnents or seciirities variously called 
tra nche s , shares, participations, classes or contracts 

25 have cash flows vrtiich cure determined by dividing emd 

distributing the cash flow of an xinderlying collection 
or pool of mortgages on a monthly basis according to 
pre- specified rules. The present value of a tranche 
can be estimated on the basis of the eaqpected monthly 

30 cash flows over the remaining term of the tranche, and 
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an estimate of the present value of a tranche can be 
represented as a xnulti -dimensional integral whose 
dimension is the nuznber of payment periods of the 
tra nch e. For a typical instrument with a 30 -year term 
5 and with monthly payments, this dimension is 360. 

Usxially, such a high- dimensional integral cam be 
evaluated only approximately, by numerical integration. 
This involves the generation of points in the domain of 
integration, evaluating or "saii5)ling" the integrand at 

10 the generated points, and combining the resulting 

integrand values, e.g., by averaging. Hell known for 
niamerical integration in securities trading is the so- 
called Monte Carlo method in which points in the domain 
of integration are generated at random. 

15 With integrands arising in financial securities 

trading, the constat ional work in combining the 
sanqpled values is negligible bb compared with producing 
the integrand values. Thus, numericeil integration 
methods in securities trading may be compcured based on 

20 the nxuhber of sanqples required for obtaining a 

sufficiently accurate approximation to the integral. 

Summary Of The Invention 

A preferred method for estimating the value of 

25 a financial security involves numerical, integration 

unlike Monte Carlo integration in that an integrand is 
sampled at deterministic points having a low- 
discrepancy property. As compared with the Monte Carlo 
m e tho d, significant advantages are realized with 

30 respect to speed, accuracy, and dependability. 

Brief Description Of The Drawing 

Pig* 1 is a schematic of a programmed conputer 
syst^ in accordance with a preferred embodiment of the 
35 invention. 
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Fig. 2 is a graphic representation of performamce 
data obtained* in computer trial nins with an exemplary 
embodiment of the invention as compared with two Monte 
Courlo computations. 
5 Further included is an i^endix with two conputer 

cQgorithms in "C source language, respectively for 
cos^ting Sobol points and Hal ton points. For a 
description of C, see B.W. Kemighan et al. , The 
Proaramr rtTTifj T^n guaqe C . Prentice- Hall, 1978. 

10 

Detail^ Description Of Preferred Embodiments 

Fig. 1 shows a stored-program computer 11 con- 
nected to input means 12, e.g., a keyboard, for enter- 
ing financial securities data, ouid connected to output 

15 means 13, e.g., a visual display device, for displaying 
an estimated value of the financial security. The 
conputer 11 includes a working memory M, a low- 
discrepcmcy deterministic point generator P, an inte- 
grand evaluator E, and an integrand - value combiner C. 

20 In estimating the value of a mult i -dimensional 

integral in financial securities trading, a 
multivariate integrand is sampled at points 
corresponding to a low-discrepancy deterministic 
sequence of points in the multivariate unit cube as 

25 defined below. If the multivariate iinit cube is also 
the domain of integration, the points of the low- 
discrepancy deterministic sequence can be used as 
sample points directly. In the case of a more general 
domain of integration, sample points correspond to 

30 points of a low- discrepancy deterministic sequence in 
the multivariate unit cube via a suitable 
transformation or mapping. 

When a sufficiently large number of. sampled values 
has been computed, an approximation of the integral is 

35 obtained by suitably combining the computed values, 
e.g., by averaging or weighted averaging. 
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In the d- dimensional iinit cube D = [0,1]**, a low- 
discrepancy deterministic sequence Z|, of 

points in D Ccui be characterized as follows: 

For a point t = [tj, — , tj in define 
5 [0, t) = [0, ti) X ... X [0, ta), 

where [0, tj) denotes a left -closed, right -open 
interval, and denote with Xto,o(») the characteristic 
or indicator function of [0, t) . For points Zj, 
Zn in D, define 
10 R.(t; Zi, z^) = (Ii=i'^ X|o,ot2k))/n - tit^.^.t^, 

cmd define the discrepancy of Zj, , as the L.- 

norm of the ftuiction R^l.; Zj, , z.) , i.e., 

|Rd(*; Zi» Zn)|«. = supigD|R«(t; Zi, z^) I . 

The sequence z,, Zj, ... is said to be a low- 
15 discrepancy deterministic sequence provided 

|r^(.; Zi, z„)i.. = O((log n)7n) . 

Low- discrepancy deterministic sequences are 
described in the published literature; see, e.g., 
H. Niederreiter, "Random Nixaiber Generation and Quasi- 
20 Monte Carlo Mfethods", CBMS-NSF . 63, SIAM, Philadelphia, 
1992 cmd S. Pckskov, "Average Ccise Complexity of 
Multivariate Integration for Smooth Functions", ^jn^^ rTiai 
of Comple xity . Vol. 9 (199J3), pp. 291-312. Well-known 
exanples of low- discrepancy deterministic sequences are 
25 the so- Celled Hammersley points, Hal ton points, Sobol 
points, and hyperbolic- cross points. 

For illustiration, in the case of Sobol points in a 
single dimension (d«l) , a constructive definition may 
be given as follows: Choose a primitive polynomial 
30 P(x) = x» + ajx^* + ... + a^,x + 1 

(whose coefficients aj are either 0 or 1) and define 
so-called direction ntunbers Vj, i > n by the following 
recurrence formula: 

Vi = aiVi,, e a^v^ e . . . e a^iVi^^., © © (v^/2°) , 
35 where © denotes a bit-by-bit "exclusive or" operation. 
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The initicd. numbers = can be 

chosen freely provided is odd and 0 < mi < 2' for 
i = 1, 2, , n. 

Using the direction numbers so defined, now 
5 define the one -dimensional Sobol sequence X|, ... 
by 

= bjVj © bjVj © • . . © bwv^, k 0 
where k « E^^rioiklb^i ±q ^he binary representation 
of k. 

10 For higher dimensions (d > 1) , the 

first d primitive polynomials Pi, Pj, P^ are 

used. If {x^^}k-i^ denotes the one -dimensional Sobol 
sequence generated by the polynomial P|, the d- 
dimensional Sobol points are defined as « (3c^\ x^, 
15 .... x^^). 

While this definition can be implemented as a 
coaster algorithm directly, faster techniques are 
known irtiich produce these points in a "shuffled" or 
permuted sequence. In particulcur, this c^lies to the 
20 con^ter algorithm given in the Appendix. 

For specificity in the following, a preferred 
stibodiment of the invention is described as applied to 
a collateralized mortgage obligation known as CMO FN, 
89-23. This has thirty-year matvirity and consists of 
25 the following tranches: 

PAC tranches 23-A, 23-B, 23-C, 23-D, 23-E 
supporting tranches 23 -G, 23 -H, 23 -J 
residual tranche 23 -R 
accrual tranche 23 -Z 
30 The monthly cash flow is divided and distributed 

according to pre -specified rules which are included in 
a fortnal prospectus- Some of the basic rules may be 
stated as follows: 

First from the monthly cash flow, the coupon is 
35 paid to the tranches. The remaining amoiint, caJ.led 
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Principal Distribution Amount , is used for repayment of 
the principal. Prior to a fixed date in the future, 
the Principal Amount will be allocated sequentially to 
the tranches 23-G, 23-H, 23-J and 23-Z. After that 
5 date, the Principal Distribution Amount will be 

allocated sequentially to the tranches 23 -A, 23 -B, 23- 
C, 23 "D and 23 -B according to a planned schedule. Any 
excess amount of the Principal Distribution Amount over 
the planned sch e dule will be eQlocated sequentieaiy to 
10 the tranches 23 -G, 23 -H, 23-J and 23-Z. A distribution 
of principal of the tranche 23 -R will be made only 
after all other tranches have been retired. 

In deriving an estimate for the present value of 
the security at the time of issue, the following 
15 notation is used below: 

C - the monthly payment on the underlying 
mortgage pool; 

ifc - the projected interest rate in month k, 
k = 1, 2, . 360? 
20 - the percentage of mortgages prepaying in 

month k; 

a^et^^i - the re m a inin g annuity after month k. 
A r emain i ng annuity cl^ can be expressed as3 to 
a^ = 1 + Vo + ... + Vo*^^ k = 1, 2, 360, 
25 with Vo = l/d+i^,), where io is the current monthly 
interest rate. Thus, after k months , the remaining 
amount of principal boinrowed is C-a^. 

It is assumed that the interest rate i^ Ccui be 
esqpressed as 
30 ik « Ko eDq>(f,)i^i, 

where exp(.) denotes exponentiation and where 
fw — / ^360 are independent, normally distributed 
random variables with mean 0 and variance a, anH ^ ±q 
a given consteuat. For the present example, a = 0.0004 
35 is chosen. 
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It is assumed fiirther that w,^ as a function of i^ 
can be computed as 

1^ ° K| -f K3 curctcm (K^ik + K4) , 
irtiere K^, K,, K^, K4 are given constcuats* 
5 Under these assumptions, the cash flow in month k, 

k = 1, 2, . 360 is 

C(l-wi (f |) ) . . . ({i, . . . , Ui) ) (1-Wk{f 1, . . w W 

+ l# • • • » fk)^36a*+l) » 

where 

^0 ^((i#--w$k) = Kj + arctan(K3Ko^ioe3^({i+. • .+(J+K4) . 

This cash flow is distributed according to the 
rules of FN, 89-23. Then, the cash flow for each of 
the tranches is multiplied by the discount factor 
Vi(fi) . . .Vk{{i, . - . ,$k) I with 
15 Vj(ei,...,fj) « 1/(1 + iyioexp(fi+...+{j)) , j«l,2, ...,360, 
to find the present va^ue for month .k. Summing of the 
present values for every month gives the present value 
PVt, for each tranche T. 

"Hie expected or estimated value » 
20 B(PVt> « B(PVT(ei,...,f36o)), 

upon a change of variables is represented by 

B(PVr) = PVyCyiCxi) , . - . ,y3«>(x3eo) - .dJCj^j, 

where 

Xi = (2ira)-**J^3fi exp(-tV(2<r))dt. 
25 Thus, for each tranche T, a 360-variate integrand 

has to be integrated over the 3 60 -dimensional unit 
cube, . 

After geneirating a point 

(X|, *3C2, ^^36o) 

30 of a low- discrepancy deterministic sequence in the unit 
cube, the point 

(Yu Yi, y36o) 
is computed by finding the value of the inverse normal 
cumulative distribution fxmction at each x^. Then, for 
35 each tranche T, the function value 
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PVx(yi/ Yt» 7360) 

Is ccsmputed. These are the function values used in 
numerical integration. 

Fig. 2 shows results from trial nms for 
5 OK) FN, 89-23 with a preferred method using Sobol 
points genes^ted by the corresponding computer 
algorithm given in the Appendix, cis conpared with Monte 
Carlo integration. Two Monte Carlo cosqcnitations were 
caanried out, with different "seeds" or steurting values 
10 of a congruent ial pseudo- random number generator known 
as RMI2; see W. Press et al Numerical Recipes in 
Cambridge University Press, 1992. It is apparent that 
the preferred method reaches a steady value more 
rapidly. In this and other trial runs, with typical 
15 requirements of precision and confidence, a speed-up by 
a factor of 3 to 5 was realized as compared with Monte 
Carlo integration. Much greater speed-up Ccui be 
expected when higher precision or/and higher levels of 
confidence eure sought. 
20 In a further trial run with CMO FN, 89-23, 

instead of Sobol points, Hal ton points were used as 
generated by the corresponding computer cd-gorithm given 
in the i^endix. It is felt that Sobol points may be 
preferred over Halton points for integrals of high 
dimension. However, this preference may not apply in 
the case of lower-dimensioncQ integrals, e.g., with 
dimension up to 5 or so. 

A coaqputatlon as described may be terminated after 
a predetermined niunber of f miction evaluations. 
Alternatively, e.g., after every function evaluation or 
after a predetermined incremental number of function 
evalxxations, a current approximation may be compared 
with one or several preceding approximations, for 
te rmin ation once a suitable condition depending on the 
difference between approximations is met. Such ter- 
mination criteria may be called automatic. Automatic 
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termination is particxilaxly reliable where a sequence 
of approximations settles down smoothly; see, e.g., 
the curve in Fig. 2 corresponding to Sobol points. 
Advantageously in computing function values, a 
5 cluster or network of multiple parallel processors or 
workstations can be used. This may involve a master or 
host processor providing points of a low- discrepancy 
sequence to slave processors and . combining function 
values returned by the slave processors into an 

10 approximate value for the integral. Thus, the crnnpu- 

tation can be speeded up in proportion to the number of 
processors used. 

Advantageous further, in combination with 
nimerical integration as described above, is the use of 

15 error reduction techniques analogous to variance reduc- 
tion in Monte Carlo integration as described, e.g., by 
M. Kalos et al., Monte Ce^^o Mgtno^g. John Wiley & 
Sons, 1986. This may involve a change of variables 
or/and variation reduction, for exanple. 
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^veid soboKiot n) 



^? i'J*^ function aebol for genezatlng Sobol pointa. It XBtuxna the 

S;^*L?:^^"»?} VOiBt, (The point (0,0 0) ia .kimed!) 

rrto point is l^pUcitly zetnrned through the ax^ x. The func^^ 
Mbeeq frea Biaarical Recipes, 1992 was used es m faaslB. but 
;^«aged eignifi^ly to s cc esttds t e the psrallel distributed s;«noach. 
The fimctica eobol can gen e rate Sobol points aUppina an iniSol 
part -of the eequenoe. The oanatant mZDns. mmm Bm»«ical Reeipea, 



ia extended to 360 which required adding nore initialiaina d^ato 
^ arrays ip (the imitiYe polynonlalB): ^^"«ig oata to 

■^Jtheir degrsM) , and iv (the initial direction nuirtiers) . me 
polynoBial x is used to generate the first coordinate of Sobol points. 

******** * . . ■ 

«. • . • • A 1 1 X « X X X XX — 

•include "nrutil.h** 
tdef ine mXBTV 30 
idef ine MAIDIM 360 

i?wf' ^ ^* di»msioo of the' points */ 

extern double •x; /* This returns istOicitly the n-th Scbol point in x */ 

void sobol (int n) 

unsigned long i,ia,ipp; 
static double fac; 

s^c unsigned long in^ixDSODlili^i) ,*iu(lsatBXT4^i] ; 
static utts&gned long ndeg|inXDiit«.lU(0,HAXaZT,l,2,3.3,4,4, 

5,5,5*5, 5v5« 

6,6,6,6,6,6, 

'^'^'7^2. 7.7,7,7,7,7,7,7,7,7,7.7,7,7, 
' " ^ „ ^ ^ ;.;r8,8,B,8,8,8,8,8, 8,8,8, 8,8,8, 

9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9, 
9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9, 

»*^» 11 •H* 11. 11- 11- n 11 51 11 

); 

/• The tint cea«Siii.t« ia met to mZBIT, but It i. nsver uaad •/ 
■tatic unaigned long ipIIJWa>lM*lJ.{0,l«XBlT,0.1,X.2,l,«, 



SUBSnrUTE sheet (rule 26) 



wo 97/07475 



PCTAJS95/103d3 



11 



/♦tog 5V 2,13,7, 14»11, 4, 

/*dog 6*/ 1,16,13,23,19,25, 

/♦deg 7*/ 1,32, 4, 8, 7, 56, 14, 28,19, 50, 21, 43, 31, 62, 37, 41, 55, 59, 

/♦deg 6*/ 14,56,21,22,38,47,49,50,52,67,70,84,97,103,115,122, 

/*&bg 9*/ 8,13 ,16 ,22 ,25,44 ,47 ,52 ,55 ,59 ,62 ,67 ,74 ,81 ,82 ,87 ,91 ,94, 
103,104,109,123,124,137,138,143,145,152,157,167,173,176,181,182,185, 
191,194,199,218,320,227,339,230,234,236,241,244,253, 

/nteg 10*/ 4,13,19,33,50,55,64,69,98,107,115,131,137,134,140,145,152, 

158,161,171,181,194,199,303,308,337,242,251,353,365,366,274,283,289,395. 
301,316,319,334,346,353,361,367,383,395,398,400,412,419,433,436,438,433, 
446 , 454 , 4 57 , 472 , 493 , 505 , 508 , 

/*degll*/ 3,11,21,22,35,49,50,56,61,70,74,79,04,88,103,104,112,115,117,122,134, 
137,146,148,157,158,162,164,168,173,185,186,191,193,199,313,214,220,227, 
236,242,251,256,259,365,266,376,393,304,310,316,319,333,338,334,339,341, 
345,346,363,367,373,375,376,381,385,388,393,409,415,416,431,438,431,434. 
439,446,451,453,457,458,471.475,470,484,493,494,499,503.517,518,534,537, 
555 , 560 , 565 , 569 , 578 , 580 , 587 , 589 , 590 , 601 , 607 , 611 , 614 , 617 , 618 , 635 , 638 , 635 . 
641,647,654,659,663.673.675,683,684,689,695,696,713,719,734,733,734,740, 
747 , 749 , 753 , 755 , 763 , 770 , 783 , 784 , 787 , 789 , 793 , 796 , 803 , 8 05 , 810 , 815 , 824 , 829 , 
- 830,832,841,847,849,861,871,878,889,892,901,908,920,933,943,949,950,954. 
961 , 968 , 971 , 973 , 979 , 983 , 986 , 998 , 1001 , 1010 , 1012 , 

/*degl2V 41,53,61,63,76,104,117,131,143,145,157,167,171,176,181,194,317,336, 
239 , 262 , 283 , 386 , 307 ) ; 

Static unsigned long iv(mZDXH*lsnZBXT+lI»{ 
0, /* 40 eleaents in a rov */ 
l,l.l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l*l.l.l,l#l,l,lrl»lrl. 

1, 1,1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,1, 1,1,1, 1,1, 1,1, 1,1, 1,1, l,l,l,l,l»l.l*l,l,lrl.lr 1*1*1. 
1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1,1, 1,1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1* 1#1.1»1»1*1*1»1»1#1«1. 
1,1,1,1, 1,1, 1,1, 1,1, 1,1, 1,1,1, 1,1, 1,1,1, 1,1, 1,1,1,1*1*1*1*1*1#1«1»1*1*1»1*1«1«1. 
1,1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1,1, 1,1,1, 1,1,1, 1,1,1, l*l»l*l*l.l*l*lrl*lrlrl*l*l*lt 
1,1,1, 1, 1, 1, 1, 1, 1, 1,1, 1,1, 1,1,1, 1,1,1, 1,1, 1,1, 1,1, 1,1, 1,1. 1*1, l»l»lt 1*1*1*1*1,1* 
.1,1, 1,1, 1.1, 1,1, 1,1. 1,1, 1,1, 1,1,1, 1,1,1,1,1,1,1*1, 1,1*1, 1*1*1*1*1*1*1*1*1*1*1*1* 
1.1.1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1*1*1*1*1*1*1*1*1*1*1*1*1*1*1* 

1,1,1,1,1,1,1,1.1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1*1*1*1*1*1*1*1*1*1*1*1*1* /*end 1 

*/ 

1,3,1,3.1,3,1,3,3,1,3,1,3,1*3,1,1,3,1,3,1,3,1*3*3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3, 

1,3,1,3,1,3,1,3.3,1,3,1,3.1,3,1,1*3,1,3,1,3,1,3,3,1,3,1,3,1,3.1,1,3.1,3,1,3,1,3, 

1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3.1,3,1.3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3, 

1,3,1,3.1.3.1.3,3,1,3,1*3,1,3,1,1,3,1,3,1,3,1,3,3,1,3*1*3,1,3*1,1,3,1,3,1,3,1,3. 

1,3.1,3.1.3.1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1.3,1,1,3,1*3,1,3,1,3, 

1,3,1,3,1,3,1,3,3,1.3.1,3,1,3,1*1,3,1,3,1,3,1,3*3,1,3,1*3*1,3,1,1*3,1*3,1,3,1,3, 

1,3,1,3,1,3,1.3,3,1,3,1,3,1,3,1,1,3,1,3,1.3.1,3,3,1,3*1,3,1,3,1,1,3,1,3,1,3,1,3, 

1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3*1*3*3,1,3,1,3,1,3,1,1.3,1,3,1,3.1,3, 

1,3,1,3,1,3,1,3,3,1,3,1*3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3*1*3.1.3.1,3, /*end 3 

*/ 

1,5,7,7,5,1,3,3.7.5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5*1,3,3.7.5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1.3,3,7,5,5,7,7,1,3*3,7,5,1,1,5,3*3,1,7,5,1,3,3.7.5,1,1,5,7,7,5,1,3.3, 

1.5,7,7,5.1,3,3,7,5,5,7,7.1,3,3,7,5*1,1*5*3*3*1*7,5,1,3,3.7,5,1,1,5.7,7,5,1,3,3, 

1.5,7,7,5,1,3,3,7,5,5,7.7,1,3,3,7,5,1,1,5,3*3*1*7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1*5*7,7,5,1,3,3,7.5,5,7,7,1,3,3,7,5,1,1,5,3,3*1*7,5,1,3,3,7,5,1,1*5*7,7,5*1,3*3, 

1*5,7,7,5,1,3,3,7,5.5.7.7,1,3,3,7,5.1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7.5,1,3.3,7,5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1.3,3,7,5,5,7,7,1,3*3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3* 

1,5,7.7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1*1,5,3,3*1,7,5,1,3,3.7,5,1,1,5,7,7,5,1.3,3, /*end 3 

•/ 

1.15,11,5,3,1,7,9,13,11,1,3.7,9,5,13,13,11.3.15,5,3.15.7,9,13.9.1.11.7.5,15,1.15,11,5,3,1,7, 
9. 

1.15,11,5.3,1,7.9.13.11,1,3,7,9,5,13,13,11,3,15,5,3,15,7,9,13,9,1,11,7,5,15,1,15,11,5,3,1,7, 
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1.15.11.5.3.1.7.9.U,ll,1.3.7.9.5.U.13.U.3.15.5.3.15.7.9.13.9.i.lX.7.5.15.HS.«.5.3.1.7. 
1.15.11.5.3.1.7.9.13.11.1.3.7.9.S.U.13.11.3.1S.5.3.15.7.9.13.9.1.11.7.5.X5.1.15.11.5.3.1.7. 
1.15.11.5.3.1.7.9.13.11,X.3.7.9.5,13,13.IX.3.1S,5.3.15.7.9.t3.9.1.1X.7.5.15.1.15.U.S.3.1.7. 
X.X5.1X,5.3.X.7.9.13,XX.X.3.7.9,5.X3.X3.Xa.3.X5.5.3.X5,7.9.13.9.1.XX.7.5.X3.X.XS.XX.5.3.X,7. 
X.15.XX.5.3.X.7.9.13.1X.X,3.7.9.5.X3.13.U,3.X5.5.3,X5.7.9.X3.9.1.XX.7.5.X5.X.X5.XX.5.3.X.7. 
X.X5.Xl,5.3.X.7.9.13.ia.X.3.7.9.5.X3.13.XX.3.15.5.3.X5.7.9.X3.9.X.X1.7.S.X5.1.XS.XX.5,3.X.7. 

3\'"x5?i:3if9f'"'''''"'"'"'"'"'"-"-«-''"-"'l'«.».3.3X.XX.5.«.27.19.2X.5.X.X7.1 
3:"l^9:3if9!'"''-''"-«'»-"-"'"-"'»-^-".".l.«.".3.3X.XX.5.«.37.X9.«.5.X.17.X 

3\'S:S.1!:S:l?:S:S:"'"'""'"'"'"'''''"'"-"'"'".»'«.».3.2x.i3.3x.25.9.«9.3 
3:S:S:;J:S:l?:J2:S:"'*'''*"'"'"'"-'-»-"'»'«'»'."'«.«^».3.2i,i3.3x.25.9.49.3 
3:S:S:5f:S:f?:iJ:S:"'"-"'"'"'"-»'''»'"'»'«.".".».a3.29.3.ax.x3.3x,a5.9.«9.3 
3':S:S:Jf;LM?;JJ:S:«'"'"'"'"'''-"'-'-»."'3»'«3.a7.x7.x5,a3.«.3.ax.x3.3x.2 

3\'S:«:;x:S.1?:JI.1J;«'"'"'"'"'^'-"'-'-*'"'39.«.a7.X7.X5.23.a9.3.ai,X3.3X,a5.9.49.3 
3:S:S:ll:SJ?:J2:S:"'"-"'"'"'"'»'''''"'39.«3.a7.X7.X5.33.29.3.ax.X3.3X.25.9.«.3 

i:S:S:ll:S:S;n:f5:"'"'"'"'"'"-"'''^'"'"«'".".«.33.a9.3.ax.x3.3x.a5.9.«,3 

3:S:«:Jl'S:|?:;J:S;«'»-»'"."'37.33.7.5,XX.39.63.27.X7.X5.23.29.3.2X.X3.3X.a5.9.49.^ 

05;|:g:SI:S!5f^ra!i9%V?o^^ 
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1,85, 67»«9, 125.25, 109,43, 89^69, X13»47«55, 97, 3»37»83a03, 37, 13, 33^15, 41,79^7. 29, 119, 75, 73,1 
05,7,59,65,21,3,113,61,89,45,107, 

1,85,67,49,125,25,109,43,89,69,113,47,55,97,3,37,83,103.27,13,33,115,41,79,17,29,119,75,73,1 
05,7,59,65,21,3,113,61,89,45.107, 

1, 85, 67, 49, 125, 25, 109, 43, 89, 69, lU, 47,55, 97, 3, 37,83, 103, 27, 13, 33, 115, 41,79, 17, 29, 119, 75, 73,1 
05,7,59,65,21,3,113,61,89,45,107, 

1,85,67,49,125,25,109,43,89,69,113,47,55,97,3,37,83,103,27,13,33,115,41,79,17,29,119,75,73,1 
05,7,59,65,21,3,113,61,89,45,107, /*eild 7 */ 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39, 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39, 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149, 15,213, 135, 253,191, 155, 175, 63,47,7,23, 39, 

1,355.79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39, 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101.255,29,203,65,195,177,105,17,47,169,109 
, 149, 15,2U, 135,253. 191, 155, 175, 63. 47, 7,23, 39, 

1.355,79,147,141,89.173,43,9,25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,353.191,155,175,63,47,7,33,39, 

1,255,79,147,141,89,173,43,9.25.115,97,19.97.197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,253.191,155,175,63,47,7,23.39, 

1.255.79.147.141.89,173,43,9.25.115.97.19,97.197.101.255.29.203,65,195,177,105,17,47.169.109 
,149, 15, 2U. 135. 253. 191, 155. 175,63, 47. 7,23. 39, 

1.255,79,147.141.89.173,43.9,25.115.97.19,97.197,101,255,29,203.65,195,177.105.17,47,169.109 
.149.15.213,135,253,191.155.175,63.47.7.23.39. /•end 8 •/ 

h?Ii!^hi^!i^^'^^'"^'^'"5'"3'*"'^^'59,353,329, 463, 385.111.475,451.263. 19, 249, 275, 
369.393,167,333,473,469,101,21,451,229,247,297,403,497.61,361, 

h?Iit"il^?iJ^'^^'"^'^^*"5,103.411,233,59,353, 329, 463. 385.111,475. 451,263, 19. 249. 275, 
369.393,167.333.473.469.101,21,451,229,247.297.403.497.61,361. 

J^Li^^il^^ ' • r 225 . 235 . 103 . 411 , 233 , 59 . 353 , 329 , 463 , 385 , 111 , 475 , 451 . 263 , 19 . 24 9 , 275 . 

369,393. 167,333.473. 469, 101, 2X.451.229. 247. 297. 403. 497, 61,361. 

i^;1*!^^f:i'^l2^^'"^'"^'*^''**^^'*"'"^'5^'^^'^^^'*"'"5,lll,475,451,263, 19. 249,275, 
369.393.1^7,333,473,469.101,21,451,229,247,397,403,497.61.361. 

385.111,475,451, 263,19,249,275, 
?^2i2'hi^i"^'*^^'*^^'"^'^^'*51,229,247,297,403,497,61,361, 

h?3At*?iJ'?4i'^''^:"^'^'^^'"^'*"'^3'59, 353,329.463,385, 111,475,451.263. 19. 349, 375, 
369.393,167,333,473,469,101.31.451.339.247,297.403.497,61,361. 

h2'Z;i^?;J^!4i'7i?^ii"'"'*"5.l03.411.233, 59.353, 329.463, 385.111.475. 451.263.19,249.375, 
369.393.167.333.473,469.101.31.451.339.347.397.403.497.61.361, 

h2^«*?^^!;i'^d*hi"'"^'*^*'"^'*"'"^'"'™'"''*"'"S.111.475.451,363,19, 349,375. 
369,393,167,333,473,469,101.31.451.329,247.297.403.497.61.361, 

h2^«J*hi^!;i'^i?^Hi®^'"5*"5'"^'*^'"^'"'^53, 329,463. 385.111.475. 451.263.19.249,275. 

/^?l2'!'^?^'"^'*^^'**''*®^*^l'*51.229.247,297.403,497.61.361. 

/•end 9 */ 

:^^*2''®''^^'^^^'"^'^^'''33.451.447.177.57,471,979, 197,251, 

;I^^;I2^;if^^!^2^f;'*!'"^'"''^"'"^'^''''"'"''"^'«57.647.581.505.833. 139, 147.203, 
,337.829,989.375.365.131,215.733.451.447.177.57.471.979.197.251. 

:*£I:"S:^*?*?^'^"'"^'^'™'*51.447.177.57,471.979, 197,251, 

:*2I:*S''^'^^'^*^*"^''^'"'3^»*51.447,177.57.471.979. 197.251. 

;n^;2^«2^;3*?^S^!;?*!'"!;*^'^"'"^'^^'*'"'"*''"'«".647.581.505,833,139,14^ 
:^!2:*S:'*?*'^'^^^'"^'*"'^5»*"'**''>^'"'57,471,979,197,251, 

;n;2^;i2^1?2:"'''*''^"^"^'^^'"'''3^'*'«"^»"5'»", 657,647,581.505,833,139.147. 203,81 
;^iZ4®S;*?!:r!i!®**"^'*^'''"'*51.447.177.57,471,979. 197,251, ' 

;^«2^;i2'i4PhPV*^^"'"^'^"'"^'^^'''«^'"''»",657.647.581.505.833. 139, 147, 203,81 
:^2Z:®5?:^*''"5'^^^'"^'^^5.733,451,447,177,57,471, 979, 197.251. 

;33;32^«-2^i''^!:"^'***'"^'^^^'^"'^57,39. 761, 169, 983, 657, 647, 581. 505.833. 139. 147. 203, 81 
,337,829.989,375,365.131,215.733,451.447,177.57,471,979.197,251. • «r * .iw.ei 
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l.rri, 731,1013, 759, 835, 949,113,939, 615, 157,39, 761, 169, 983,657, 647, 561,505. 833, 139, 147, 203,81 

,337,839,969,375,365,131,215,733,451,447,177,57,471,979,197,251, 

/*ettd 10 •/ 

1,1265,023,727,267,833,471,1601,1341,913,1725,2021,1905,375,893,1599,415,605,819,975,915,171 

5,1223,1367,663,629,525,469,981,1667,1587,1251,451,461,721,483,1209,1457,415,1435, 

1,1285,823,727,267,833,471,1601,1341,913,1725,2021,1905,375,893,1599,415,605,619,975,915,171 

5,1223,1367,663,629,525,469,961,1667,1587,1251,451,461,721,483,1209,1457,415,1435, 

1,1285, 823, 727, 267»833,471,1601, 1341,913, 1725,2021, 1905, 375,693, 1599,415,605,619,975,915,171 

5,1223,1367,663,629,535,469,981,1667,1567,1251,451,461,721,463,1209,1457,415,1435, 

1,1265,623,727,267,833,471,1601,1341,913,1725,2021,1905,375,893,1599,415,605,819,975,915,171 

5,1223,1367,663,629,525,469,981,1667,1567,1^1,451,481,721,463,1209,1457,415,1435, 

1,1265,823,727,267,833,471,1601,1341,913,1725,2021,1905,375,893,1599,415,605,819,975,915,171 

5,1223,1367,663,629,525,469,981,1667,1567,1251,451,461,721,483,1209,1457,415,1435, 

1,1285,823,727,267,833,471,1601,1341,913,1725,2021,1905,375,093,1599,415,605,619,975,915,171 

5,1223,1367,663,629,525.469,981,1667,1587,1251,451,461,721,483,1209,1457,415,1435, 

1,1265,823,727,267,633,471,1601,1341,913,1725,2021,1905,375,893,1599,415,605,819,975,915,171 

5,1223,1367,663,629,525,469,981,1667,1587,1251,451,481,721,483,1209,1457,415,1435, 

1,1285,823,737,267,833,471,1601,1341,913,1725,2021,1905,375,093,1599,415,605,819,975,915,171 

5 , 1223 , 1367 , 663 , 629 , 525, 469 , 981 , 1667 , 1587 , 1251 , 451 , 481 , 721 , 463 , 1209 , 1457 , 4 15 , 1435 , 

1 , 1265 , 823 , 727 , 267 , 833 , 471 ,1601 , 1341 , 913 , 1725 , 2 021 , 1905 , 3 75 , 693 , 1599 , 415 , 605 , 619 , 975 , 915 , 171 

5,1223,1367,663,629,525,469,981,1667.1587,1251,451,481,721,483,1209,1457,415,1435, 
/*eBd 11 */ 

1,3855,4091,987,1839,4033,2515,579,3863,977,3463,2909.3379,1349.3739,347,387,2881,2821,1873, 
1959 , 1939 , 2389 , 3251 , 1149 , 243 , 3609 , 1131 , 1701 , 143 , 1339 , 3497 , 2499 , 1571 , 9 83 , 4031 , 1635 , 3317 , 1163 , 



for (j-ndeglkli-ljjc.ilAXBrr.-j'M.) { 
ipp*ip(k) ; 
i<=iu(j-iBdeg[k)n)c),- 
i - (i » ndegCk]) ; 
for (l-»de9(X)-l;l»l;l--) ( 
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if Cipp & i> 1 iutj-1) M; 
»• 1; 

iuti)m.i; 

) * 

fae>l. 0/(11. « MMEBm; 



) 

else 

/* Chedc if the (n-1) -th nuHber was ganexated in the previous call 

to mtsbol. It BOt, update in and ix */ 

i£(inl-n-l) { 
unaignM long gny; 

/* Set ix to 0 */ 

for (k»l;k<»IlIZB(d«MAZDXH) ;lc-M-) ixPcJ-O; 
iniA-1; 

gray^in^ (in»l) ; /* rind gray code of in */ 
for (j-lrjc^oaBXTpj-M-) ( 

if (gxsy&l) { /• Only digits whiai are 1 are used */ 

iai»<j-l)*ISKXDXM; 

for (l&«l;lc<>IMIBCd,MAZDIll);k-M>) ixM ivfimfk] ; 
^ gx«y»«l; 

laain; /* Calculate the next -vector in the 

for (j»l;j<-aiAXBrr;j-M'> I /* rind the rigfatnost aero hit •/ 

if (Mia 1)) break; 

in »• 1; 

if (j > MAXBTD nrerzorCllAXBXT too bmII in sobseq**) ; 
iB»(j-l)*linZDXM; 

for (k«lrk<«ZlIIH(d*llAXDZM);k+^) ( 
ixPtl iv(iflM^k); 
x{k-ll»ix[k)*fae; 

) 

MMDZT 
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This is tfaa fmiction tmlton foe gftnexati&g Hal ton points. 

It retuxBS tbs n-th d-diaasionsl Hal ton point. The point is iaplicitly 

returned through the array x. The laat two digits of n-l in base p(j] 

are kq^ in ql(j] and <|2[jl. Wien both digits beoone p[j] -1, the 

radical inverse function is oosvuted again. TSiat my the aecunulntion 

of round*off error is acvoided. In practice, there are not any upper bounds 

on the values of d and n. 



extern int d; /* actual diMOsion of the points */ 

extern int *4l,*q3f /* qlti) is the last digit of n-l in base p[jl, q2[jl is 

the digit before the laat one */ 
extern double *x; /* This vill contain the Balton point */ 
extern int *p; /* the first d priM m^ers V 
extern int *p_l; /* fixet d priae nuebers minus l */ 

extern double *incn,*lncr2; /• inenfj] is l/plj) and l/(p(j] *p[jl > */ 

double find_fi(int p, int n) t /* See belov */ 

void balton < int n) 

double a; 
int j»nn; 

static int ins.n; /* ibe default value of insjn is 0 */ 

/* Che^ if the (n-l) -tb nnitor was generated in the previoos call 
to halton. If not» update qi» 42, and x */ 

if |ins_ntan-l) 

ins_n>n-l; 
^ for <j-0; j<df j++> 

4l[j)-ins.n%p[j]; 
«2CjJ-(ins_n/pljJ »%p[j) ; 
^ x(j]^ind_fi(p(j],ins_n); 

ins 

) ~ 
for<j.P;j<d;j++) 

^ i«CqlIj)^_ltjI) 

/* Xt is easy to update when the last digit is less than p(j) -l */ 

qiiji+^; 

^ xfj)"XCj)-i-incn(j]; 
^ else if C«J2(J)<p_i[jJ) 

/* Ibis is the case when the last digit is pH] -1 and the digit 
before the last one is leas than pf jl -1 */ 

qltj]-0; 

^ x(j]sx{j]-i'incrl[j]-i'incr2[j]-1.0; 
else 
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/* This ia the case lAen the last digit ia p[j) -1 and the digit 
before the last one is also p(j)-X */ 

n»-»/«PWl*pCJ))f 

ai«m%ptj); 

im-nn/pbJ; 

if (im) xtj]-(a4>£iad fi(p(j),i]n))*iaex3Cj]*lncnCj]; 
^ else X [ j] ^nncra t jl^^iacn t j] ; 

) ^. 
/* 

find.fidnt int n> 



double £ind_fl(iiit p, int n) 

ixkt p2«nii; 
double 8,fi»isicr; 

ittcr»l.O/p; 

iBMl/pa; 



while (nn > 0) 

ti*m (o&%p)*a; 



^ s*ainer; 

/* me two largest rfpm i ei i ts of f i are added later to avoid possible 

loss of predsioa */ 
fi«C(a/p%p)^ fi>/ p2; 
fi-f(n%p)/ (double) p; 
^ return fi; 



SUBSTITUTE SHEET (RULE 26) 



wo 97/07475 PCT/US95/l()3d3 

Claims 

1 1. A compater method in financial securities trading, 

2 for producing an approximate value for cm estimated 

3 value of a financial security, con^rising: 

4 providing the ccmputer with financial security 

5 data; 

6 producing the approximate value by nuxaerically 

7 integrating a multivariate integrand whose multi- 

8 dimensional integral over a domain of integration 

9 represents the estimated vcaue, wherein nxmterical 

10 integration comprises ; 

11 evaluating the integrand at points in the 

12 domain of integration corresponding to points of a 

13 low- discrepancy deterministic sequence, and 

14 combining the integrand values to produce the 

15 approximate value; and 

16 producing the approximate value as an output 

17 for inspection or/and fuirther processing. 

1 2. The computer method of claim l, ^^rein combining 

2 the integraxid values comprises averaging. 

1 3. The computer method of claim 1, wherein the nuxhber 

2 of integrand values is predetermined. 

1 4. The computer method of claim 1, wherein the ntunber 

2 of integrand values is determined automatically . 

1 5. The computer method of claim 1, further craqprising 

2 allocating integrand evaluations among a plurality 

3 of processors. 

1 6. The computer method of claim 1, further comprising 

2 application of an error reduction technique. 
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1 7. The computer method of claim 6, wherein error 

2 reduction comprises a change of variables. 

1 8. The computer method of claim 6, wherein error 

2 reduction comprises variation reduction. 

1 9. The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence conprises Sobol 

3 points. 

1 10* The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence comprises Halton 

3 points. 

1 11. The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence comprises 

3 Hammersley points. 

1 12. The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence comprises 

3 hyperbolic -cross points. 

1 13. The computer method of claim 1, wherein financial 

2 securities data comprise derivative instnnnent data. 

1 14. The computer method of claim 1, further comprising 

2 tising the approximate value in offering the security 

3 for sale. 

1 15. The computer method of claim 1, further conqprising 

2 using the approximate value in deciding whether to 

3 buy, sell or hold the security. 

1 16. A coaster system for financial securities trading, 

2 for producing an approximate value for an estimated 

3 value of a financial security, comprising: 
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4 means for providing t:he cox^niter with financial 

5 security data; 

6 tnp?"« for producing the approximate value by 

7 numerically integrating a multivariate integrand 

8 whose multi-dimensional integral over a domain of 

9 integration represents the estimated value, 

10 cooqprising: 

11 Tnopwc for evalxiating the integrand at points in 

12 the domain of integration corresponding to points of 

13 a low* discrepancy deterministic sequence, and 

14 means for combining the integrand values to 

15 produce the approximate value; and 

16 means for producing the approximate value as an 

17 output for inspection or/and further processing. 

1 17. The computer system of claim 16, wherein the means 

2 for combining the integrand values comprises 

3 averaging means. 

1 18. The computer systoi of claim 16, wherein the nimtber 

2 of integrand values is predetermined . 

1 19. The computer system of claim 16, further conprising 

2 means for automatically determining the number of 

3 integrand values . 

1 20. The computer system of claim 16, further comprising 

2 means for collocating integrand evaluations among a 

3 plurality of processors* 

1 21. The computer syst^ of claim 16, fiirther conqprising 

2 means for applying an error reduction technique. 

1 22. The cdt^ter system of claim 21, wherein the means 

2 for applying an error reduction technique comprises 

3 means for a change of variables. 
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1 23 . The ccsmputer system of claim 21, whezrein the t^p 

2 for applying an error reduction technique comprises 

3 means for variation reduction. 

1 24. The computer syst^ of claim 16, wherein the. means 

2 for producing a low-discrepcmcy deterministic 

3 secpience comprises means for producing Sobol points 

1 25. The computer system of claim 16, wherein the tno^inff 

2 for producing a low-discrepancy deterministic 

3 sequence comprises means for producing Halton 

4 points . 

1 26. The computer system of claim 16, wherein the means 

2 for producing a low*discrex>ancy deterministic 

3 sequence comprises means for producing Haimnersley 

4 points . 

1 27. The computer syst^ of claim 16, wherein the means 

2 for producing a low-discrepancy deterministic 

3 sequence comprises means for producing hyperbolic - 

4 cross points. 
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